Nonlinear driven diffusive systems with dissipation: fluctuating hydrodynamics 
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We consider a general class of nonlinear diffusive models with bulk dissipation and boundary driv- 
ing, and derive its hydrodynamic description in the large size limit. Both the average macroscopic 
behavior and the fluctuating properties of the hydrodynamic fields are obtained from the micro- 
scopic dynamics. This analysis yields a fluctuating balance equation for the local energy density at 
the mesoscopic level, characterized by two terms: (i) a diffusive term, with a current that fluctuates 
around its average behavior given by nonlinear Fourier's law, and (ii) a dissipation term which is a 
general function of the local energy density. The quasi-elasticity of microscopic dynamics, required 
in order to have a nontrivial competition between diffusion and dissipation in the macroscopic limit, 
implies a noiseless dissipation term in the balance equation, so dissipation fluctuations are enslaved 
to those of the density field. The microscopic complexity is thus condensed in just three trans- 
port coefficients, the diffusivity, the mobility and a new dissipation coefficient, which are explicitly 
calculated within a local equilibrium approximation. Interestingly, the diffusivity and mobility co- 
efficients obey an Einstein relation despite the fully nonequilibrium character of the problem. The 
general theory here presented is applied to a particular albeit broad family of systems, the sim- 
plest nonlinear dissipative variant of the so-called KMP model for heat transport. The theoretical 
predictions are compared to extensive numerical simulations, and an excellent agreement is found. 



I. INTRODUCTION 

In many different contexts like biological physics, ac- 
tive matter, electronics, combustion, granular media, 
population dynamics or chemical reactions, to name just 
a few, the dynamics of the system of interest is char- 
acterized at the mesoscopic level by reaction-diffusion 
equations [1-9] . This type of equations often arises from 
the competition between diffusion and dissipation mech- 
anisms, which typically drives the system out of equilib- 
rium. The dissipative character of the dynamics implies 
a continuous loss of energy to the environment, so an 
steady energy input is needed in order to maintain the 
system in a stationary state, which is usually attained by 
a boundary injection mechanism. In addition, most sys- 
tems in this class are strongly nonlinear, with transport 
coefficients which depend on the local energy density. 

The physics of nonlinear driven dissipative systems is 
poorly understood as a result of several factors. On one 
hand this family of systems is intrinsically out of equilib- 
rium: there is no equivalent to the Gibbs distribution of 
equilibrium systems to describe the statistics of micro- 
scopic configurations. Furthermore, there are currents of 
mass or energy induced by local gradients, which in turn 
are controlled by the dissipation and not by boundary 
conditions, as in standard, conservative nonequilibrium 
systems. In addition, microscopic dynamics in driven 
dissipative media is typically irreversible, which leads to 
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complications in their statistical description. 

Despite these difficulties, there have been recent ad- 
vances in nonequilibrium physics which are opening new 
avenues of research [TUl E] , offering tools to understand 
the macroscopic behavior of driven dissipative media 
starting from their microscopic dynamics [T2J [T3] . The 
key idea which has triggered these developments has been 
the realization of the essential role played by macro- 
scopic fluctuations, both in equilibrium and away from 
it. In fact, the study of fluctuation statistics of macro- 
scopic observables provides an alternative way to obtain 
thermodynamic potentials, complementary to the usual 
ensemble description. This observation, valid both in 
equilibrium [14] and nonequilibrium |101 111] , is however 
most relevant in the latter case, where no general bottom- 
up approach connecting microscopic dynamics to macro- 
scopic properties has been found yet. The large deviation 
function (LDF) controlling the statistics of these fluctua- 
tions plays in nonequilibrium statistical mechanics a role 
similar to the equilibrium free energy |15[ I16j . and the 
computation of LDFs in different nonequilibrium systems 
has thus become one of the main objectives of nonequi- 
librium statistical physics in recent years, triggering an 
enormous research effort which has led to some remark- 
able results [Ifi l im|T5H2T l . Th e calculation of LDFs from 
first principles (i.e. from microscopic dynamics) is typi- 
cally a daunting task, which has been accomplished only 
for a handful of models, in most cases simple interacting 
lattice gases [TDJ[TT]. There is however an alternative the- 
oretical framework, the Macroscopic Fluctuation Theory 
(MFT) of Bertini, De Sole, Gabrielli, Jona-Lasinio and 
Landim [10] , which studies dynamic fluctuations in diffu- 
sive media at a mesoscopic level, and offers explicit pre- 
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dictions for LDFs arbitrarily far from equilibrium. The 
starting point for this theoretical scheme is the fluctuat- 
ing hydro dynamic equations describing at the mesoscopic 
level the evolution of the system of interest. From these 
equations, using a path integral formulation, one may ob- 
tain the probability of paths in mesoscopic phase space 
associated to a given fluctuation of a macroscopic observ- 
able, and from those probabilities a variational problem 
for the optimal path responsible of a given fluctuation 
and the associated LDF is derived [H2 EJ [El Q2H2]. 

The above advances have been mostly restricted to 
conservative nonequilibrium systems, where dissipation 
is absent and nonequilibrium conditions are solely in- 
duced via boundary gradients or external fields (see how- 
ever refs. [THUS])- In order to extend the ideas of Macro- 
scopic Fluctuation Theory to nonlinear driven dissipative 
systems, it is of utmost importance to develop minimal 
models for this broad class of systems which, while cap- 
turing their essential ingredients (namely nonlinear diffu- 
sion, bulk dissipation and boundary driving), are simple 
enough to be amenable to both analytical calculations 
and extensive computer simulations. In particular, as 
the starting point of MFT is the mesoscopic hydrody- 
namic description of the system at hand, it is essential to 
obtain the fluctuating hydrodynamic equations and the 
associated transport coefficients. 

In this work, we analyze the hydrodynamic behavior of 
a general class of systems whose main ingredients are dif- 
fusion, dissipation and boundary driving, which together 
with nonlinear behavior are the fingerprints of many re- 
alistic driven dissipative media. In our family of models 
there is one particle at each site of a d-dimensional lattice, 
and the state of each particle is completely characterized 
by its "energy" . The dynamics is stochastic and pro- 
ceeds via collisions between nearest-neighbor particles, 
at a rate which is a general function of the pair energy. 
In a collision, a certain fraction of the pair energy is dissi- 
pated, and the remaining energy is randomly distributed 
between the two particles. In addition, the system may 
be coupled to boundary thermal baths. 

In the large system size limit, both continuous space 
and time variables can be introduced, as well as the "hy- 
drodynamic" fields: energy density p(x, t), current j(x, t) 
and dissipation d(x,t). The time evolution of the energy 
density follows a fluctuating balance equation of the form 

d t p(x,t) = -d x j(x,t) + d(x,t) . 

Interestingly, the microscopic dynamics must be quasi- 
elastic in order to ensure that both diffusion and dissipa- 
tion take place over the same time scale in the continuum 
limit. Using a local equilibrium approximation, the cur- 
rent and dissipation fields can be expressed as functions 
of the local energy density. In particular, the fluctuat- 
ing current can be written as j = —D(p)d x p + £, where 
the first term is nothing but Fourier's law with a diffu- 
sivity D(p), and the second term £ is a noise perturba- 
tion, white and gaussian. The current noise amplitude is 
a(p)/L (i.e., the noise strength scales as L~ x / 2 ), where 



L is the system size and a(p) is often referred to as the 
mobility in the literature. These gaussian fluctuations 
are expected to emerge for most situations in the appro- 
priate mesoscopic limit as a result of a central limit theo- 
rem. Microscopic interactions can be highly complicated, 
but the ensuing fluctuations of the slow hydrodynamic 
fields result from the sum of an enormous amount of ran- 
dom events at the microscale which give rise to Gaussian 
statistics at the mesoscale. In the present case, a proof 
of the gaussian character of the noise can be given, due 
to the simplicity of the class of models considered. On 
the other hand, the dissipation field can be written as 
d = —i/R(p), where v is a macroscopic dissipation coeffi- 
cient which can be related to the inelasticity of the under- 
lying microscopic dynamics, and R(p) is a new transport 
coefficient. The dissipation field has no intrinsic noise, 
so its observed fluctuations are enslaved to those of the 
density. This stems from the subdominant role of the 
noise affecting the dissipative term: its strength scales as 
L~ 3 / 2 as a consequence of the quasi-elasticity of the mi- 
croscopic dynamics, so it is negligible against the current 
noise in the mesoscopic limit. 

The class of models introduced in this paper repre- 
sents at a coarse-grained level the physics of many dis- 
sipative systems of technological as well as theoretical 
interest. In particular, when the collision rate is con- 
stant (i.e. independent of the pair energy), the Kipnis- 
Marchioro-Presutti (KMP) model [25] for heat conduc- 
tion is recovered in the conservative limit. The KMP 
model plays a main role in nonequilibrium statistical 
physics as a benchmark to test theoretical advances [T71 - 
12*2*] . For the dissipative case, different generalizations of 
the KMP model have been recently proposed [El [23] . 
Our general class of models contains the essential ingre- 
dients characterizing most dissipative media, namely: (i) 
nonlinear diffusive dynamics, (ii) bulk dissipation, and 
(iii) boundary injection. Moreover, it can be regarded as 
a toy model for dense granular media: particles cannot 
freely move but may collide with their nearest neighbors, 
losing a fraction of the pair energy and exchanging the 
rest thereof randomly. The inelasticity parameter can be 
thus considered as the analogue to the restitution coeffi- 
cient in granular systems [BJ. 

In a forthcoming paper [24] we will use the fluctuat- 
ing hydrodynamic picture that emerges from this work 
as starting point to analyze the large deviation function 
of the dissipated energy for this general class of nonlin- 
ear driven dissipative systems. In order to do so, we will 
extend the tools of Macroscopic Fluctuation Theory to 
this broad class of systems, and test our results in ad- 
vanced Monte Carlo simulations capable of probing the 
rare events associated to the tails of the LDFs of interest. 

The remainder of paper is organized as follows. In Sec- 
tion II we define the general class of nonlinear driven dis- 
sipative models mentioned above and discuss some limits 
of interest. Section III is devoted to the derivation of the 
hydrodynamic evolution equation for this general family 
of models, which is of reaction-diffusion type. In Section 
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IV we study the fluctuating corrections to the hydrody- 
namic equation, deriving in this way the fluctuating hy- 
drodynamics for our family of models, from which a full 
characterization of their large deviations statistics can be 
obtained via Macroscopic Fluctuation Theory jTDl [51] . It 
is shown that both the current and dissipation noises are 
white and gaussian, but the dissipation noise is subdom- 
inant against the current noise in the large system size 
limit. In Section V we study in detail a particular family 
of models, for which the collision rate is a power of the 
local energy. We solve the full hydrodynamic problem, 
and compare our results with extensive numerical sim- 
ulations, finding excellent agreement. Furthermore, this 
agreement extends also to the transport coefficients and 
the fluctuations of the hydrodynamic fields. A summary 
of the main results of the paper, together with a physi- 



cal discussion thereof, is given in sec. VI The appendices 
deal with some technical details that, for the sake of clar- 
ity, we have preferred to omit in the main text. 



II. A GENERAL CLASS OF NONLINEAR 
DRIVEN DISSIPATIVE MODELS 

Let us consider a general class of models whose main 
ingredients are diffusion, dissipation and boundary driv- 
ing. For the sake of simplicity, we will present them for 
the one-dimensional (ID) case, but the extension to ar- 
bitrary dimension is straightforward. The system is thus 
defined in a ID lattice with N sites. A configuration at 
a given time step p is given by p = {pi tP }, I = 1, . . . , N, 
where pi iP > is the energy of the l-th site at time 
p. Thus, the total energy of the system at this time is 
E p = ^2fi—xPi,p- The dynamics is stochastic and proceeds 
as follows. In an elementary step, a nearest neighbor pair 
of sites (1,1 + 1) "collides" with probability 

/(£*,p) 



j i,p — pi,p + pi+i,p> 



(2.1) 



where / is a given function of the pair energy S; iP , and L 
is the number of possible pairs. Clearly L ~ TV, but the 
particular relation depends on the boundary conditions 
(e.g., L — N + 1 for open boundary conditions while 
L = N for the periodic case). Once a pair is chosen, 
a certain fraction of its energy, namely (1 — a)S; jP , is 
dissipated to the environment. The remaining energy 
aS/ iP is then randomly redistributed between both sites, 

Pi, P +i = z p aT,i tP , pi +ltP+1 = (1 - z p )aT,i }P , (2.2) 

with z p an homogeneously distributed random number in 
the interval [0, 1]. This dynamics defines the evolution of 
all bulk pairs, / = 1, . . . , N — 1. In addition, and depend- 
ing on the boundary conditions imposed, boundary sites 
might interact with thermal baths at both ends, possibly 
at different temperatures Tl (left) and Tr (right). In this 
case the dynamics is 

Pi.p+i = z P a(e ltP + e L ), Pn, p +i = z p a(e NtP + e R ), 

(2.3) 
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FIG. 1. The model is defined on lattice sites, each one charac- 
terized by an energy pi . The dynamics is stochastic and pro- 
ceeds via random collisions between nearest neighbors where 
part of the pair energy is dissipated to the environment and 
the rest is randomly redistributed within the pair. Such dy- 
namics mimics at the mesoscopic level the evolution of a wide 
class of systems characterized by a nontrivial competition be- 
tween diffusion and dissipation. 



when the first (last) site interacts with its neighbor- 
ing thermal reservoir. Here e Vl v — L,R, is an en- 
ergy randomly drawn at each step from the canonical 
distribution at temperature T„, i.e. with probability 
prob(e u ) = T~ x exp(— e v /T v ) (our unit of temperature 
is fixed by making ks — 1), see Fig. [I] We may consider 
instead an isolated system with periodic boundary condi- 
tions, such that L = N and Eqs. ( |2.1[ ) and (2.2) remain 
valid for I = (I = N) with the substitutions po p = pn, p 

{PN+l,p = Pl.p)- 

The simplest dynamics corresponds to /(Ej p ) = 1 in 
Eq. (2.1 1. In this case all (nearest neighbor) pairs collide 
with equal probability P; jP — 1 1 



L , independently of their 
energy. This choice (together with a = 1 above) corre- 
sponds to the Kipnis-Marchioro-Presutti (KMP) model 
of heat conduction 22 , which can be considered as a 
coarse-grained description of a large class of ID diffusive 
systems of technological and theoretical interest. It also 
plays a main role in non-equilibrium statistical mechanics 
as a benchmark to prove rigorous results and test theo- 
retical advances. For instance, it is one of the very few in- 
stances where Fourier's law can be rigorously proved |22) . 
In addition, the KMP model has been used to investigate 
the validity of the additivity principle for current fluctu- 
ations [19] and the Gallavotti-Cohen fluctuation theorem 
[T7] and its generalization in refs. [201 HI] ■ Another sim- 
ple, but physically relevant, choice is /(£;. p ) — E;, p , so 
that Pi tP ~ Y±i iP /(2E p ) for a large system. A variant 
of this model has been recently used to study compact 
wave propagation in microscopic nonlinear diffusion [25] . 
In general, as the total exit rate from configuration p - 
the denominator of Eq. (2.1)- is a sum of L terms, we 
expect that P; ;P oc L —1 , so it is convenient to write 



PiAp) 



m p (L) : 



with 



1 



tt p (L)^ z ]Tj-(£ ; , p ) 
i=i 



(2.4) 

so £lp(L) remains finite as L — > oo. Throughout this sec- 
tion, we will analyze the stochastic process generated by 
the dynamics ( 2. 1 )-( 2.2 ) supplemented with the appro- 



priate boundary conditions. 

Note that energy is conserved in the dynamics only 
for a = 1, while it is continuously dissipated for any 
< a < 1. In this way, the parameter a can be consid- 
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ered as the equivalent of a restitution coefficient, using 
the language of granular materials. Thus, in an isolated 
system (without boundary driving) the energy would de- 
crease monotonically in time. However, if energy is in- 
jected, for instance via coupling to boundary thermal 
baths as described above, a steady state will be eventu- 
ally reached where energy injection and dissipation bal- 
ance each other. In the conservative case a = 1, the 
transition rates introduced above verify detailed balance 
with the "microcanonical" distribution, i.e. the distribu- 
tion giving equal probabilities to all the microstates com- 
patible with the total energy of the system E. This im- 
plies that a closed system would eventually approach the 
microcanonical equilibrium distribution for any choice of 
the function / appearing in the transition rates. There- 
fore, a large subsystem, comprising a large number of 
sites, would be described by the canonical distribution 



with its temperature imposed by the rest of the system, 
which acts as a thermal reservoir. Thus, if we take a 
part of the system containing a large number of sites, 
it is expected on a physical basis that this region would 
approach the "local equilibrium" state compatible with 
the average values of the relevant macroscopic quantities, 
if the latter follow a very slow dynamics as compared 
to the relaxation of the system to the local equilibrium 
state. Afterwards, over the much slower "hydrodynamic" 
scale, the relevant macroscopic quantities evolve accord- 
ing to the macroscopic equations and reach a steady state 
compatible with the boundary conditions imposed on the 
system. 



III. HYDRODYNAMIC DESCRIPTION 



The dynamics defined in Eq. (2.2) can be rewritten as 
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(3.1) 
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for each realization k (trajectory) of the stochastic pro- equation (3.1) can be rewritten in the following way, 



cess. Here Sij is Kronecker's delta, 



» 



is a random 



number homogeneously distributed in [0, 1] which con- 
trols the local energy exchange, and y p <E [1, L] is an 
independent random integer which selects the colliding 
pair. The (conditional) probability for y p being equal to 
I, provided that the system is in a certain configuration 
p, is given by Pt tP (p), see Eqs.(2.1) and (2.4), 



y P ,0p 



m p (L) 



(3.2) 



The energy current involved in a collision of a pair 
(I, I + 1), denoted here by ji tP , is just the net energy trav- 
eling to the right in an unit time once dissipation takes 
place. Mathematically, again for each trajectory k of the 
system, 
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(3.3) 



where we used Eq. (2.2) in the second equality. Equiva- 



lent^, we define the dissipated energy at site I at time p 
over each trajectory as 



H,p 



■(!-«)/©( 



Vp 



i +S v 1 p"\i- 



(3.4) 



when the colliding pair is either (I, I + 1) or (I — 1, 1). No- 
tice that, as opposed to the current term, the dissipation 
is defined for each site I instead that for each pair (I, l+l). 
This definition simplifies the analysis below. Making use 
of the current and dissipation definitions, the evolution 



Pl,p+x Pi,p - Jl-i,p h, P + a i, P ' y 6 -^) 

that is a discrete reaction-diffusion equation, with a dif- 
fusive contribution ji—i tP — ji,p and a sink term d^ p . The 
latter is proportional to 1 — a, thus vanishing in the elas- 
tic case a = 1. 

The hydrodynamic (or average) evolution equation is 
obtained by summing over all the possible trajectories 
of the system k, i.e. if A is a certain physical property 
of the system, (A) = Y.IZi^ / n t, where A^ is the 
value of A in the fcth realization of the stochastic pro- 
cess, and the number of trajectories Nt — > oo. This 
procedure is equivalent to averaging over all the possible 
sequences of the pair of independent random numbers 
(z p , y p ) which determine the microscopic dynamics. This 
said, and in order not to clutter our formulas, we will 
drop the superindex k in the remainder of the paper, as 
all expressions where the pair of random numbers (z p , y p ) 
appears are also valid for each trajectory of the system. 
Averaging in this way Eq. (3.5) we arrive at 



(pl p+ i - pi. p ) = (ji-i, P - 3l,p) + (di, P ) 



(3.6) 



Let us consider first the average value of the energy cur- 
rent. From Eq. (3.3), 



/■ v _ a / {Pi,p - Pl+l,p)f(El,p) 
[Jl ' p> 2L\ fip(L) 



(3.7) 



where we have taken into account Eq. (3.2) and (z p ) — 
1/2. On the other hand, the average value of the dissi- 
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pation is 



(di, P ) = L " ',1 



/(S i|P )+/(S t - 1 , p ) 

n p (L) 



(3.8) 



In the large system size L — > oo, we will be interested 
in the density, current and dissipation fields, that are 
smooth functions of the continuous spatial variable 



/ 

L 



Ax 



Xl + \ - Xi 



1 

V 



(3.9) 



with x G [—1/2, 1/2]. In this continuum limit, the aver- 
age density (pi, p ) will be replaced by p av (x,t), 



(Pl,p) Pav(x,t), 



(3.10) 



where t is a continuous time that will be introduced later 
on. On the other hand, because of the discrete derivative 



inside the average in (3.7), (ji p ) is expected to scale as 



(jl,p) -> L 2 T P j^v(x,t), 



r p = lim (^(L)). 



(3.11a) 



(3.11b) 



We have introduced the term t p above for the sake of 



convenience, because (f2~ 1 (L)) has a finite limit as L — > 
oo. Since for each configuration of the system Q p is its 
total exit rate, t p is a sort of microscopic time scale, that 
depends on the choice of the collision rate function /(£). 
Finally, the dissipative term (d; jP ) in Eq. (3.8 1 scales 
as (1 — a)L _1 . From this discussion, the diffusive term 



in the evolution equation (3.6) for the average density, 
should scale as L~ 3 , 



i-i, P - jLp) 



i-u 



Ji, 



-3, 



p^xjav^X 7 £) 



(3.12) 



and this diffusive term should be neglected in the ther- 
modynamic limit, unless the "microscopic" dissipation 
parameter scales as 1 — a oc L~ 2 . In fact, this is the 
correct scaling for the microscopic "inelasticity" 1 — a 
in the large system size limit, being the only one which 
guarantees that both diffusion and dissipation interplay 
on the same time scale at the mesoscopic level 12] . With 
this scaling, 



(di, P ) 



L 3 T n d : 



p u av {x, t) , 



1 



V 

2L 2 '' 



(3.13) 



where v is a "mesoscopic" dissipation coefficient which 
remains finite in the large system size limit L — > oo. 

In order to be more concrete, we have to evaluate the 
averages on the rhs of Eqs. (3.7) and (3.8). This can 



be done within a local equilibrium approximation (see 
appendix [A]) , which we expect to be valid in the ther- 
modynamic limit [TD1|2B]. For large system sizes, a clear 
time scale separation appears between the scale (of the 
order of several times t p ) in which the system approaches 
locally the equilibrium distribution (i.e. a local prod- 
uct Gibbs measure with local temperature (pi lP )) and the 



much longer time scale over which the average density 
field evolves following the hydrodynamic equation. Un- 
der this local equilibrium approximation, the average cur- 
rent verifies 



(jl, P ) L 2 T p D({pi iP )) 



(pi+i, P - Pi, P ) 



with the diffusivity D(p) given by 



D(P) 



6 



dr r 7 f(pr 2 )e 



(3.14) 



(3.15) 



Equation (3.14) is thus consistent with the scaling in Eq. 



(3.11a). Furthermore, 

3&v{Xi tj 



D(p av )d x p 

av ' 



(3.16) 



i.e. Fourier's law holds with diffusivity D(p av ). On the 
other hand, the average dissipation, Eq. ( |3.12 ), can be 
written within the local equilibrium approximation as 



with 



(di, P ) ~ -L T p vR((pi yP }), 



R{p) = p / drr b f(pr z )e- 
Jo 



(3.17) 



(3.18) 



Again, Eq. (3.17) is consistent with the previously as- 
sumed scaling (3.13), with 



d av (x,t) = -vR(p av ) 



(3.19) 



Interestingly, this new transport coefficient R{p), associ- 
ated to the dissipation, can be related to the diffusivity. 
By differentiating Eq. (3.18) with respect to p after a 



change of variables z = r^/p, it is found that 
1 dR{p) R{p) 



D{p) 



6 dp 



3p 



(3.20) 



In fact, given the diffusivity, this equation can be consid- 
ered as a first order differential equation for R(p). The 
solution thereof, with an appropriate boundary condition 
(normally R(p = 0) = 0), is equivalent to calculate the 
integral in Eq. (3.18). 



Now, by introducing a macroscopic time t, such that 
its time increment at step p is given by 



p-i 



p-i 



At p = L-% 



(3.21) 



we can rewrite the evolution equation for the average 
density (3.6) as 



dtp a v(x, t) = -d x j av (x, t) + d av (x, t), 



(3.22) 



where the average current j av (x,t) verifies Fourier's law 
with diffusivity D(p av ), Eqs. ( |3.15| )-(3.16l, and the aver- 
age dissipation d av ( x, t) is given in terms of the density 
field by Eqs. (0^-(|3T9l. 
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The boundary conditions for eq (3.22 1 depend on the IV 



physical situation of interest. For instance, we may con- 
sider that the system is kept in contact with two thermal 
reservoirs at x — ±1/2, at the same temperature T. In 
that case, the system eventually reaches a steady state 
in the long time limit, for which the injection of energy 
through the boundaries and the dissipation balance each 
other. The stationary average (macroscopic) solution of 
([3~22l verifies 



JavO) + vR{Pw{x)) = 0, j av (x) = -D(p av (x))p' av (x), 

(3.23) 

where the prime indicates spatial derivative. The first 
equation in (3.23) follows from (3.22), while the second 



one is Fourier's law. Equivalently, a closed second-order 
equation for p may be written, 



fix 



[D(p av )p' av ] = uR(p av ), 



(3.24) 



with the boundary conditions p av (±l/2) = T. We may 
now introduce an auxiliary field y that will be helpful 
later on, 



such that 



y = R(j>), 



d{x,t) = -i/y(x,t), 



(3.25a) 



(3.25b) 



Equations (3.23) and (3.24) can be also rewritten for the 
variable y, 

jL v ( x ) + vy&v(x) = 0, j av (x) = -D(y av (x))y' av (x), 

(3.26) 

with 



D(y) = 



±y i D{p)= (*my i D 



dp 



\ dp J 



(p) . (3.27) 



Thus, D acts as an "effective" diffusivity, i.e. the factor 
multiplying the spatial gradient when writing Fourier's 
equation in terms of the new variable y. Taking into 
account Eq. (13. 20}, 



D{y) 



1 1 d\np(y) 
6 3 d In y 



(3.28) 



Interestingly, the above equation shows that D is con- 
stant, independent of y, whenever y = R{p) depends al- 
gebraically on p. This observation will be very useful in 
Section |Vj where we study in detail a particular family 
of models. Equation (3.26) can also be summarized in a 



second order differential equation for y av , 

D(y, v )y' av ] ' = vy av , y«,(±i/2) = R(T), (3.29) 

the solution of which gives the average density field in 
the situation at hand. When D is constant, Equation 



(3.29) is linear in y and can be readily integrated. 



FLUCTUATIONS OF THE CURRENT AND 
DISSIPATION FIELDS 



The local energy p^ p obeys the balance equation (3.5), 
with a diffusive term, given by the discrete spatial deriva- 
tive of the microscopic current, and a dissipative term. 
Taking averages over all possible realizations of the dy- 
namics, the hydrodynamic equation (3.22) has been ob- 



tained, with a similar structure but now in terms of the 
the continuum hydrodynamic fields (density, current and 
dissipation). The main objective of this section is to an- 
alyze the local fluctuations of the current and the dissi- 
pation fields away from their averages or, in other words, 
the properties of their respective noise terms. 

We will address the problem of characterizing the cur- 
rent and dissipation noise terms in the limit of large sys- 
tem size L » 1, the same one in which Eq. (3.22) has 



been shown to be valid over the hydrodynamic scale de- 
fined in Eq. ( |3.21 1. The idea is to split the microscopic 



current j; iP into a "main" term j/ iP , whose average co- 



incides with (ji tP ), see Eq. (3.7), and a noise term £; p . 
Similarly, we will have g?; iP = di p + r\i p , with the main 
term d^ p verifying (d^ p ) = (fi; iP ), see Eq. (3.8). Thus, 
r\i iP is the dissipation noise. We will investigate the prop- 
erties of both noise terms and find that, in the large sys- 
tem size limit: (i) both noises are white and Gaussian, 
and (ii) the dissipation noise is subdominant agai nst th e 
current noise if 1 — a = 0(L~ 2 ), as given by Eq. (3.13). 
Therefore, the fluctuations of the dissipation field are en- 
slaved to those of the density field, and the current noise 
is the main source for the fluctuations in the system at 
the mesoscopic level. 



A. Fluctuating current. Average value and noise 
properties 

Let us consider the fluctuations of the microscopic cur- 
rent j'z jP . As we have already stated, the idea is to split 
the current into a "main" and a "noise" term, 



3i,p — 3i,p + €l,p 



(4.1) 



such that (ji, p ) = (ji,p)- It is clear from Eq. (3.7) that 
the following choice 



Ji, P = 



a (pi, p - Pi+i,p) /(Sj, P ) 



2L 



(4.2) 



guarantees the previous condition on the average. To 
stress out the difference between jt tP and ji tP , notice that 
while the former is exactly zero unless pair (1,1 + 1) ac- 
tually collides at time p, see Eq. (3.3), the latter may 



take a nontrivial value even in the absence of such col- 
lision; however their averages coincide [27]. Eq. (4.2) is 
Fourier's law at the microscopic level, i.e. ji jP is propor- 
tional to the local density gradient. We want to study 
now the properties of the noise £; iP = j/ jP — ji, p , namely 
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its average and its correlation function Since 
(ji,p) = (jl,p), it is c l ear that 



(4.3) 



For the noise correlation function, it is straightforward 
to show that 



<6,p6',p'} = (ji.pji'.p') - (ji.pji'.p')- 



(4.4) 



Furthermore, from the current definition Eq. (3.3), it 
is also easily proven that (Ci,p^i',p') = for p ^ p' , so 
the current noise at different times is uncorrelated. For 



equal times, p — p' , the second term on the rhs of (4.4) 



is negligible against the first one in the limit L > 1, 
because it is ~ (D(L~ 2 ), while the leading behavior of 
the first term will be shown to be ~ 0(L _1 ). Using now 
the definition of the microscopic current (13.31) , 



jl,pjl',p = a 2 [(1 - z p )pi tP - Zppi+i, p ] (4.5) 
x [(1 - z v )p v . p - z p pi>+i tP ] S Vpi iSy p M. 

Taking into account that (z p ) = 1/2, (z 2 ) = 1/3, and 

m p {L) M '' 



yp,l°Vp,l'/P 



(4.6) 



we arrive at 



{jl,p3l',p) 



of 
Si' 



(pi+i,p + Pi, P " Pl,pPi+i,p) ' 



n p (L) 



(4.7) 

Hence, the leading behavior of the current noise correla- 
tion in the L>1 limit is 

1 . 



\£i,p£i',p') ~ ySi, p 8ij>S p y, (4.8) 
where we have neglected (D(L~ 2 ) terms, and defined 



[Pi+i,p + Pi, P ~ Pi,pPi+i,p) /( s ;,p) ' 



n p (L) 



(4.9) 



Equation (4.9) cannot be expressed in terms of (pi iP ), so 
the fluctuating balance equation (3.5) is not closed at 



the microscopic level, as otherwise expected. Using here 
again the local equilibrium approximation previously em- 
ployed [TU] (see appendix |A|) , we obtain 



2 2 



drr 7 f(p av r 2 )e~ 



(4.10) 



A rigorous proof of this local equilibrium approximation 
can be found in [26], in the context of conservative in- 
teracting particle systems. The application thereof to 
the dissipative case studied here relies on the fact that 
the microscopic dynamics is quasielastic. Moreover, the 
validity of this local equilibrium approximation will be 
fully confirmed a posteriori in extensive numerical simu- 
lations, see Section [Vj 



In the large system size limit we have to take into ac- 
count Eq. (3.11a), and introduce a similar scaling for the 



fluctuating current and its noise, 

j Wp -> L- 2 T p j(x, t), Cl, P -> L- 2 t p £(x, t), (4.11) 
so that 



j(x,t)=j(x,t)+£(x,t). 



With this definition, 



(j(x,t)) = jav(M), 



(4.12) 



(4.13) 



and thus (£(x,i)) = 0. Besides, the correl ation functi on 
of the current noise is, by combining Eqs. (4.8 )-(4.11 ), 



(£(x, t)£(a/, f )) ~ ^(p &v )5(x - x')6(t - if), (4.14) 

where we have taken into account the quasielasticity of 
the microscopic dynamics, 1 — a = 0(L~ 2 ), as given 



by Eq. (3.13). We have also introduced the so-called 
mobility 



Jo 

and made use of 

Sim 



drr 7 f(pr 2 )e- r2 =2p 2 D(p) (4.15) 



Ax 



At r . 



(4.16) 



where Ax and At p were defined in Eqs. (3. 9) and (3.21), 



respectively. The current noise is hence white, its av- 
erage value vanishes and it is delta-correlated in space. 
Furthermore, it can be shown that it is gaussian, and 
we present a sketch of the proof in appendix [5] Remark- 
ably, Eq. (4.15), which relates mobility and diffusivity, 



a(p) = 2p 2 D(p), is a fluctuation-dissipation relation for 
the dissipative case. In fact, it is the same one as in the 
conservative case, because of the quasi-elasticity of the 
underlying microscopic (stochastic) dynamics 28J. 

The gaussian character of the current field fluctuations 
will allow us, in a forthcoming paper [24] . to extend the 
recently-introduced Macroscopic Fluctuation Theory [10] 
to the case of nonlinear driven dissipative systems [13] . 
This theoretical scheme allows to study the fluctuations 
(both typical and rare) of macroscopic observables arbi- 
trarily far from equilibrium, offering predictions for the 
large-deviation functions which characterize their statis- 
tics and the optimal paths in phase space responsible of 
these fluctuations. The validity of the fluctuating hydro- 
dynamic equation derived in this paper, and in particu- 
lar of the local equilibrium approximation used here, can 
be checked a posteriori by comparing the predictions for 
the large-deviation functions that arise from the macro- 
scopic fluctuation theory associated to our equation with 
the numerical results. 



B. Fluctuating dissipation. Average value and 
noise properties 



Our starting point is Eq. (3.4) for the microscopic 



dissipation d/ iP , that we split into a main term plus a 
noise term, i.e. 



di, p — d^ p + rji_ p 



1 - « Pl, P [f(E>i, P ) + /( s i-i,p)] 



L 



n p {L) 



(4.17b) 



Clearly the average of the main term verifies (di tP ) — 
(dip), see Eq. (3.8), and thus the noise average vanishes, 



(Vl,p) = 0- From Eqs. (3.4) and (4.17), we have 

m,p = —(i— a )pi,p 



'Vp, 



+ 5 y v -i-i 



+ /(Ej-ij,) 



m p {L) 



(4.17a) and 



(4.18) 



Vi&Vl' ,p' = (l-a) 2 Pi,pP;', P ' 



/(S i|P ) + /(S t -i, P ) 

lq p (l) 



(4.19) 



Now, taking averages it is clear that (rji^riii jP >) = if p ^= p', because the variables y p , y p / are independent in that 
case and the terms in brackets have zero average. Then, we restrict ourselves to the case p — p' ', for which a simple 
calculation similar to the one carried out for the current leads to 



<M',p) ~ ^^(ttTTV lf&l, P )(hi' + 8v>-i) + m-ij,M,i, + <W)] ) 



Q p {L) 



(4.20) 



after having made use of Eq. ( |4.6[ ) and neglected (D(L~ 2 ) 
terms. In order to give an explicit expression for this 
correlation, we should have to evaluate the averages in 



( 4.20 ) , for instance by making use of the local equilibrium 



approximation, as in the previous subsection. Neverthe- 
less, for the work presented here, this is not necessary 
due to the quasi-elasticity of the microscopic dynamics, 
as given by Eq. (3.13). This will allow us to show that 



the dissipation noise is subdominant against the current 
noise in the large system size limit, without having to 
calculate explicitly its amplitude. The structure of (4.20 ) 
implies that 



{Vl, P m',p') ~ L 



.(2); 



H,V 



>i,v+i 



PiP 1 
(4.21) 



where 



.(0 



= 1,2,3, are certain averages which remain 



of the order of unity in the limit L — > oo (they become 
functions of (pi, p ) in the local equilibrium approxima- 
tion). W ith th e consistent definition r)(x,t) = L 3 T]^ p /t p , 
see Eq. (3.13), we introduce the continuum limit of the 



dissipation noise and its correlation function becomes 
(r}(x, t)r){x', t')) - L- 3 v 2 k(p)5{x - x')5{t - t'), (4.22) 



where we have used Eq. (4.16), and defined n(p) as the 



continuum limit of ftJ^"above. In particular, it can 
be shown t hat th e continuum limit of the diagonal term 
Kjp in Eq. (4.21 ) is exactly equal to a(p), while each non- 

diagonal term n\ contributes a(p)/A to n(p). Equation 
(4.22) tells us that the dissipation noise (oc L~ 3 / 2 ) is 



much weaker than the current noise (oc L^ 1 / 2 ) in the 
large system size limit L> 1. Therefore, we will neglect 
it in the following, i.e. we consider that 



d(x, t) = —vR{p), 



(4.23) 



i.e. only the fluctuations affecting the current term will 
play a relevant role in the mesoscopic limit. We expect 
that this approximation becomes exact in the thermody- 
namic limit as L — ¥ oo. Then, the fluctuations of the dis- 
sipation d(x, t) are "enslaved" to those of the energy field 
p(x, t), as a consequence of the quasi-elasticity of the dy- 
namics at the microscopic level. This quasi-elastic char- 
acter of the microscopic dynamics is compatible with the 
existence of a finite dissipation at the mesoscopic level, 
as expressed by the finite value of the macroscopic dissi- 
pation coefficient v. 



V. A NONLINEAR DISSIPATIVE VERSION OF 
THE KMP MODEL OF HEAT TRANSPORT 

We now will apply the theory developed in the previous 
sections to a broad class of dissipative models. More con- 
cretely, we will restrict ourselves to the following choice 
for the collision rate function 



f(p) = 



1X8 + 3)' 



(5.1) 



i.e., f{p) oc p 13 , with j3 > —3 but otherwise arbitrary. 
We have introduced the constant 2/T(/3 + 3) [29] for the 
sake of convenience, as it simplifies the expressions of the 
transport coefficients, see below. For f3 = 0, f(p) — 1 and 
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all the pairs collide with equal probability, independently 
of their energy value. Thus, the dissipative generalization 
of the KMP model introduced in ref. [T3] is recovered. 
For (3 = 1, f(p) = p/3 and the colliding pairs are cho- 
sen with probability proportional to their energy. The 
conservative case has been recently analyzed in |25j . 

The transport coefficients for this family of models are 
easily calculated. The coefficient linked to the dissipa- 
tion, R(p), is readily obtained from Eq. (3.18), 



R(j>) 



T{fi + 3) 



P + l / drr 5+2P ( 



(5.2) 



that gives the rationale behind the choice of the propor- 
tionality constant in Eq. (5.1). The diffusivity is calcu- 



lated by substituting Eq. (5.4) into (3.20) 



D(p) = 



(3 + 3 



(5.3) 



Of course, the same result is obtained with Eq.(3.15) 



Finally, the mobility u{p) follows from the fluctuation- 
dissipation relation (4.15), which gives it in terms of the 
diffusivity, 



a(p) = 2p 2 D(p) 



(3 + 3 



,/3+2 



(5.4) 



Of course, for (3 = the values of the transport coeffi- 
cients of the dissipative version of the KMP model are 
recovered, D(p) = 1/2, a(p) = p 2 , and R(p) = p [13]. 
Interestingly, this kind of dependence with the energy 
density p appears in real systems. For instance, in gran- 
ular materials [5] the density field p may be assimilated to 
the local granular temperature. Moreover, for the hard 
sphere model, the average collision rate is proportional 
to the square root of the granular temperature. Thus, 
this granular gas case should correspond to (3 = 1/2, and 
in fact it is found that D(p) oc p 1 / 2 while the dissipative 
term goes as R(p) oc p 3 / 2 . The latter is the responsi- 
ble for the algebraic decay of the granular temperature 
(Haff 's law; oc t~ 2 for large times) observed in the homo- 
geneous case when the system is isolated [5] . 

For the class of models at hand, we thus have a fluc- 
tuating hydrodynamic equation 

d t p{x, t) = -d x j(x, t) - vR(p(x, t)), (5.5) 

where j(x,t) is the fluctuating current, 

j(x,t) = Q[p(x,t)]+Z(x,t), (5.6) 

with local average behavior given by Fourier's law, 



Q[p} = -D{p)d x p{x,t), 



(5.7) 



and v is the macroscopic dissipation coefficient defined 
in Eq. \3.\3\. The first term in the rhs of Eq. 



accounts for the diffusive spreading of the energy, and it 
is also present in the conservative case, while the second 
one gives the rate of dissipation of energy in the bulk. 




FIG. 2. Decay with time of the average energy per site in the 
isolated nonlinear dissipative KMP model. Here we plot data 
for v = 1, L = 80, and different values of (3 £ [0,2], though 
qualitatively similar behavior is obtained for other values of 
v and L. Lines are theoretical predictions, see Eq. (5.9). 



The current noise term £(cc, t) describes the random fluc- 
tuations at the mesoscopic level, as has been analyzed 
in the previous section. In the large system size limit 
L 3> 1, it is gaussian and white with 

(£(!,*)> =0, (5.8a) 



(C(as, mx\ 0) = l°(p)S(x - x'W - t'). (5.8b) 

This gaussian fluctuating field emerges in the mesoscopic 
limit as a result of the central limit theorem: although 
microscopic interactions are rather complicated, the en- 
suing fluctuations of the hydrodynamic fields, which 
evolve over a much slower time scale, result from the 
sum of an enormous amount of random events at the mi- 
croscale and give rise to gaussian statistics of C(L -1 / 2 ) 
at the mesoscale. 

Let us investigate now the average behavior for the 
family of models defined by the collision rate choice (5.1 1, 
whose transport coefficients are given by Eqs. (5.2)-(5.4). 



First, we focus on the time evolution of an isolated sys- 
tem (i.e. with periodic boundary conditions). Due to 
the dissipation, and assuming no spatial structure, the 
system enters into an homogeneous cooling state (HCS) 
characterized by a continuous energy dissipation. Ho- 
mogeneity implies that p(x,t) = p(t), and thi s en ergy 
dens ity obeys dtp = —vp^ +1 according to Eqs. (5.5) and 
(5.2). The solution is 



(3 = 



10 



where po is the initial energy density in the system (which 
is fixed) . Notice the exponential decay of energy in time 
for the linear (/? = 0) case, while for all nonlinear cases 
f3 7^ a power-law decay is obtained. This is in fact 
a generalization of the well-known Haff 's law of granu- 
lar gases, which corresponds to j3 = 1/2 in our model. 
Fig. [2] shows the time evolution of the average energy 
density for different /?'s and particular values of v and 
L as measured in Monte Carlo simulations. The agree- 



ment between numerical results and Eq. ( 5.9 1 is excellent 



in all cases, supporting the homogeneity conjecture used 
above. We do not expect the HCS to become unsta- 
ble as time proceeds because our model does not contain 
velocity variables whose correlation upon collision is at 
the origin of the well-known breakdown of HCS in d- 
dimensional (d > 1) granular gases [BJ. 

Next, we study the steady-state behavior when the sys- 
tem is coupled to boundary thermal baths at equal tem- 
perature T. Following the procedure introduced at the 
end of Section [nj it is convenient to define an auxiliary 
variable 



y : 



(5.10) 



where we have made use of Eqs. (3.25a) and \h.2\. We 



also calculate the "effective" diffusivity D(y) defined in 
Eq. (p7|, 



D 



13 + 3 



(5.11) 



where we have made use of Eq. (3.281. Thus, while the 
"tru e" diffusivity D(p) depends on p, as given by Eq. 
( |5.2[ ), the "effective" diffusivity is constant, D(y) = D. 
This allows us to calculate explicitly the average profiles 
for the density and the current, since Eq. (3.291, that 
determines them, is linear for y. The steady average 
solution in this boundary-driven case is thus 



cosh 



T 



13+1 



cosh 




(5.12a) 



sinh 



j av (x) = -Dy' &v (x) = -TP+WvD 



cosh 



(5.12b) 

The average density is readily obtained by combining 
([5~T0]) and ((5~l2l, 



p av {x) = T 



cosh | x yi -g 
cosh . 



(5.13) 



The average dissipation field is basically y av (x), since 

d av (x) = -i/y av (x). (5-14) 




FIG. 3. Average density profiles for different values of v g 
[1CT 1 , 10 3 ] and ,3 = 0, 0.5, 1, for N = 160. The bath temper- 
ature is T = 1. Lines correspond to theoretical predictions, 
and the agreement is excellent in all cases. 



Equation (5.131 reveals clearly one of the main prop- 
erties of dissipative media: the gradients are controlled 
by the dissipation and not by the boundary conditions. 
Although both ends of the system are in contact with two 
heat reservoirs at the same temperature T, the density 
is not constant throughout the system, as it is clearly 
shown in Fig. [3] the higher the macroscopic dissipation 
coefficient v gets, the larger the spatial gradient is. In the 
figure, the points correspond to Monte Carlo simulations 
of the stochastic model introduced in Sec. |TlJ while the 
lines are the analytical prediction of the hydrodynamic 
theory, Eq. (5.13). The agreement between theory and 
simulation is excellent in all cases. In the quasielastic 
limit v — > 0, the density is almost constant, as p'(x) oc v 
in that limit. This is shown in the figure by the case 
v = lO^ 1 , in which a weak spatial structure is observed in 
the scale of the figure, for all values of f3. As v increases, 
the system departs from the "quasi-homogeneous" be- 
havior. In fact, for the highest value v = 10 3 , there are 
two boundary layers near the system walls, characterized 
by a very high internal gradient, while p — > in the bulk 
of the system. 

An important question is the rate of convergence to 
the hydrodynamic description as the system size is in- 
creased. In order give an estimation thereof, we plot in 
Fig. [ijthe difference Apl(x) = pi^l(x) — p av (x) between 
the numerical values measured for the density and the 
theoretical ones for different system sizes in the range 
10 < L < 160, as a function of the spatial coordinate 
x. For the sake of concreteness, we show the data for 
v = 10 and (3 — 0.5, for which the spatial gradient is 
quite important in the system, as seen in Fig. [3j Ap^ (x) 
exhibits a nontrivial structure which however rapidly de- 
creases with the system size, being very small across the 
whole system for L = 160. In fact, even for the smallest 
size L = 10, ApL is quite small, since the relative error in 
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FIG. 4. Average excess density profiles for different values of 
L 6 [10,160], v = 10 and /3 = 0.5. Finite size corrections 
quickly decay to zero with the system size. 
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FIG. 5. Boundary thermal resistance as a function of the 
system size L for different values of v 6 [10 -3 , 10 2 ] and j3 = 
0.5, 1. In all cases the boundary thermal gap decays as L -1 , 
with an amplitude which increases with v. 



under one per cent for all x. This means that the hydro- 
dynamic description developed here is valid already for 
quite small system sizes, although the concrete threshold 
size depends on the macroscopic dissipation coefficient v. 
As v increases, larger system sizes must be considered, 
since the gradients become more important. 

A main feature of the nonlinear driven dissipative mod- 
els studied here is the presence of a nontrivial bound- 
ary thermal resistance which shows up as a gap between 
the average energy of boundary sites and the tempera- 
ture of the corresponding heat baths, A7Y = ph^,(x — 
±|) — T 7^ 0, see Figs. [3]and|4] Interestingly, this bound- 
ary thermal resistance appears only for /3 ^ 0, meaning 
that the nonlinearity is essential to develop a boundary 
thermal gap. This phenomenon is well-known to appear 



FIG. 6. Average current profiles for different values of v £ 
[10~\l0 3 ] and p = for N = 160. Notice the nontrivial 
structure due to bulk dissipation (in contrast to conserva- 
tive systems where the current field is constant across space) . 
Lines correspond to theoretical predictions, and the agree- 
ment is excellent in all cases. 



in nonlinear energy transport as a result of the scatter- 
ing of energy carriers when crossing the bath interface, 
with examples ranging from Fermi-Pasta-Ulam oscillator 
chains to hard-disks fluids, etc. [3U]. Fig. [5] plots the 
measured boundary thermal gap as a function of the sys- 
tem size for different values of v and j3 > 0, showing that 
AT^ decays as L^ 1 for large enough systems, with an 
amplitude that increases with the dissipation parameter 
v and the nonlinearity /3. 

Further insight into the hydrodynamic behavior dis- 
cussed above is given by Fig. |6| The average current pro- 
file, as measured in simulations after each collision event, 
is shown for different values of the macroscopic dissipa- 
tion coefficient v. For the sake of concreteness, and to 
make the figure appearance simpler, we have considered 
/3 = 0, though other values of j3 exhibit qualitatively sim- 
ilar behavior. The current is not constant throughout the 
system, in contrast to the situation found in conservative 
systems. For high dissipation, v ^> 1, the average cur- 
rent vanishes except in the boundary layers, inside which 
it is in fact very large. Again, the agreement between 
the theoretical expression for the current (solid line), Eq. 
(5.12b), and simulation (points) is excellent in all cases. 



The spatial dependence of the average dissipation field 
davix) is presented in Fig. [7] Note the logarithmic scale 
in the vertical axis, in order to show more clearly the 
behavior in the bulk of the system in the high dissipa- 
tion regime v ^> 1. Again the agreement between the- 
ory (solid line) and simulation (points) is excellent in 
all cases, up to the very center of the system, in which 
|rfav(£)| is really very small when v 1. 

The simplicity of the model at hand allows us to mea- 
sure directly in simulations the transport coefficients en- 
tering the fluctuating hydrodynamic theory. In particu- 
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FIG. 7. Average dissipation profiles for different values of 
v G [10"\10 3 ] and ,3 = 0, 0.5, 1, for N = 160. The bath 
temperature is T = 1. Lines correspond to theoretical predic- 
tions, and the agreement is excellent in all cases. 




FIG. 9. Comparison between the numerical and theoretical 
values of the diffusivity D(p) for j3 — 0, 0.5, 1. Different val- 
ues of v and L have been considered. The agreement between 
theory and simulation is excellent. 




FIG. 8. Dissipation transport coefficient R(p) measured in 
simulations as a function of the local average density, for dif- 
ferent values of v, j3 and L. For each value of (3, the numerical 
data collapse on the theoretical prediction p^ +1 W, L. 



lar, a strong implication of this hydrodynamic descrip- 
tion is that the transport coefficients are functions of 
the energy density p alone, for large system sizes and 
under the local equilibrium hypothesis, as expressed by 
Eqs. (13. 151), (13. 181) and (|4.15|. We have tested the valid- 



ity of this picture by comparing the numerical values of 
the transport coefficients, measured in Monte Carlo sim- 
ulations, and the theoretical predictions for the general 
class of models considered in this section, Eqs. ( |5.2[ ) and 
(5.3). Let us start by considering the numerical evalua- 
tion of R(p), the new transport coefficient associated to 



the dissipation. We have built Fig. [8] by plotting the 
local dissipation along the chain versus the local energy 
density, i.e. by combining the points of Figs. [3] and [7J 
and plotting \d av (x)\/u vs p av (x) Vx 6 [— Aside 
from the system size L — 160 considered in Figs. [3] and 
[7j we have included points obtained for smaller system 
sizes, ranging from L = 10 to L — 80, in order to see the 
convergence to the hydrodynamic behavior as the system 
size grows. For each value of /?, all the points collapse 
onto a master curve, which agrees perfectly with the the- 
oretical expression (5.2 1. This is clear signature of the 



rapid convergence to the hydrodynamic description with 
the system size, L = 10 being enough to recover the the- 
oretical prediction. 

In order to measure numerically the diffusivity, which 
is just the factor relating the local energy current with 
the local density gradient according to Fourier's law, we 
have to combine the data used in Figs. [3] and [6] to plot 
~jav/Pav (where the prime indicates the numerical spa- 
tial derivative of the measured energy profile) vs the av- 
erage energy density p av . As for the dissipation transport 
coefficient the agreement between simulation and 

theory is again very good for D{p). Notice however that 
there are some points which depart from the the theoret- 
ical prediction, but this departure is due to the smallness 
of both j av (x) and p av (x) near the center of the system, 
which gives rise to large numerical errors for the quotient 

— jav/ Pav 

Finally, we have also checked the theoretical prediction 
for the noise amplitudes, by comparing them to the nu- 
merical data obtained in Monte Carlo simulations. Fig- 
ure [lO] shows the measured amplitude of the current fluc- 
tuations, the so-called mobility a. Combination of Eqs. 



13 




p=05 




P=1 .0 



13=0.5 



p=0.0 



FIG. 10. Comparison between the numerical measurement of 
the mobility <r(p), as given by Eq. (5.151, and the theoretical 
prediction ( |5.4[ ), for /3 — 0, 0.5, 1. The agreement is excel- 
lent in all cases, even for rather small system sizes (L = 10). 
Curves for /3 = 0.5 (/3 = 1) have been shifted vertically by 0.3 
(0.6) units. The inset shows the equivalent plot in logarithmic 
scale. 



(OL Q and p~15| gives that 



O-(Pav)- 



(5.15) 



Therefore, we have plotted the lhs of the above equation 
as a function of the average density along the chain, both 
of them numerically measured in simulations. This have 
been done for the same values of /3 of the previous fig- 
ures, and for various values of v and L. In the inset, 
the same plot is shown in logarithmic scales. Straight 
lines with slope j3 + 2 are observed, in agreement with 
the algebraic behavior of a o c p@ +2 predicted by the hy- 
drodynamic theory, Eq. (5.4 1. It should be stressed that 



the excellent agreement between theory and simulation, 
and particularly the collapse of the data corresponding 
to different values of L, is a strong proof that the pre- 
dicted scaling of the current noise with the s ystem size, 
(£(x,t)£(x', If)) oc L^ 1 as given by Eq. (4.14), is correct. 
Similarly, Fig. 11 shows the measured diagonal fluc- 



tuations of the dissipation field as a function of the local 
average density. Following Eqs. (4.21 )-(4.22 ) and the 



ensuing discussion, these fluctuations are defined as 



Kdiag(Pav) = Cr(p av ). 



(5.16) 



Again the agreement between theory and simulation re- 
sults is very good, though finite-size effects, particularly 
at high densities, are more evident here that for the cur- 
rent fluctuations, see Fig. [To} This is natural since dissi- 
pation fluctuations are very small and have to be scaled 
as L 5 in the continuum limit, which strongly amplifies 



FIG. 11. Comparison between the numerical measurement 
of the diagonal part of n(p), a s gi ven by Eq 
the theoretical prediction Eq. 



(5.161, and 
(Ok, for /? = 0, 0.5, 1 and 



different values of v and L. Curves for /3 = 0.5 (/3 = 1) have 
been shifted vertically by 0.3 (0.6) units for clarity. 



the small finite-size effects present in measurements. In 
any case, Fig. [IT] confirms the scaling of the dissipa- 
tion fluctuations, i.e. that (r)(x,t)r](x' ,t')) oc L~ 3 in the 
continuum limit, as given by Eq. (4.22). This supports 



our claim that dissipation fluctuations are negligible in 
the L — > oo limit, where the only source of fluctuations 
concerns the current field. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have studied in detail the fluctuating 
hydrodynamic theory for a large class of diffusive lat- 
tice models characterized by the presence of bulk dissi- 
pation and boundary driving. The models dynamics is 
stochastic: pairs of neighboring lattice sites collide at a 
rate which is a certain function of the pair energy. A 
given part of the pair energy is dissipated to the en- 
vironment, and the remainder thereof is randomly re- 
distributed within the pair. In addition, boundary sites 
may be coupled to thermal reservoirs which on average 
inject energy into the system so eventually it reaches a 
(nonequilibrium) steady state. Thus, this class of driven 
dissipative models is a generalization of the KMP model 
for heat conduction to the nonlinear dissipative case. 

In the large system size limit, the energy density p(x, t) 
has been shown to obey a mesoscopic fluctuating balance 
equation, which contains both fluctuating current and 
dissipation fields, j(x,t) and d(x,t), respectively. The 
average current field obeys locally a non-linear Fourier's 
law, with a diffusion coefficient which depends on the 
density, while the average dissipation becomes a certain 
function of the local density. Moreover, the current noise 
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is shown to be white and gaussian, while the dissipation 
noise is subdominant because of the quasi-elasticity of 
the microscopic dynamics. Detailed expressions for the 
transport coefficients are derived in this context, using a 
local equilibrium approximation whose validity has been 
tested a posteriori. Interestingly, the amplitude o~(p) of 
the current fluctuations (or mobility) is related to the 
diffusivity D(p) by an Einstein fluctuation-dissipation re- 
lation, <j(p) — 2p 2 D(p), despite the fully nonequilibrium 
character of the problem at hand. Besides, the new trans- 
port coefficient associated to the dissipation, R{p), is not 
independent of the diffusion coefficient, but it is related 
thereto by a first order differential equation. 

In order to test in detail the emerging picture, we have 
applied the general fluctuating hydrodynamic framework 
developed in the first part of the paper to a particular 
family of models, in which the collision rate depends al- 
gebraically on the energy of the chosen pair. The rea- 
sons for this choice are twofold: (i) such power-law de- 
pendence mimics the physics of many realistic systems 
(e.g. for hard-disk fluids the collision rate is propor- 
tional to the square root of the local energy density), and 
(ii) with this choice, the transport coefficients can be ex- 
plicitly calculated as simple, also algebraic, functions of 
the energy density. We have checked the predictions of 
the fluctuating hydrodynamic theory against extensive 
Monte Carlo simulations of the proposed family of mod- 
els, and an excellent agreement is found, even for quite 
small system sizes L > 10. This agreement strongly sup- 
ports the validity of the fluctuating hydrodynamic pic- 
ture and, given the generality of the proposed family of 
models, a similar situation should be expected for any 
non-pathological choice of the collision rate function. 

The fluctuating hydrodynamic equation here derived 
is a first step in the complete description of the macro- 
scopic behavior for the broad class of nonequilibrium 
models introduced in this paper. The governing evolution 



law, summarized in just three transport coefficients, can 
now be used in conjunction with the recently-introduced 
Macroscopic Fluctuation Theory (MFT) [TU] to derive 
the large-deviation functions (LDFs) which control the 
statistics of fluctuations of the relevant macroscopic ob- 
servables, in this case the current and the dissipation. 
LDFs play in nonequilibrium systems a role akin to the 
free energy or entropy function in equilibrium systems, 
and hence are of utmost importance to the development 
of nonequilibrium statistical physics [10] . In a forthcom- 
ing paper [23] (see also [T3]) we will use the tools of MFT 
and advanced Monte Carlo simulations to study in de- 
tail the LDF of the relevant macroscopic observables for 
the general class of nonlinear driven dissipative systems 
here introduced, thus aiming at a full characterization of 
macroscopic behavior for this broad family of systems. 



Appendix A: Local equilibrium approximation 

In the hydrodynamic fluctuating theory developed in 
this paper, one of the main goals in to write down a closed 
fluctuating balance equation for the energy density at 
the mesoscopic level. At some points in the development 
of this theoretical framework, we have found that the 
average current, Eq. (3.7), the average dissipation, Eq. 



( 3.8 ) , or the amplitude of the current noise, Eq. ( 4.9 ) , are 



given in terms of certain average values with the general 
form 



g{pi ,pi Pi+i,p 



(Al) 



where g(pi tP , Pi+i.p) are different functions of the energies 
{pi, p , while S; iP and Q P (L) have been defined in 

Eq. (2.1) and (2.4), respectively. Therefore, 



' g(Pl,p, Pl+i,p)f(Pl,p + Pl+l,p) 



oo poo 



L V\-^) I JO 

pOC pOO 



dpi I dpi +1 g(pi,pi +1 )f(pi + pi +1 )P(pi,p i+1 ;p) 



n p (L) 

x/ dpi--- I dpi-x I dpt + 2--- I dp N n~ 1 (L)P(p 1 ,...,pi- 1} pi +2 ,...,pN\pi,Pi+i]p)- 

(A2) 



In order to write the equation above, we have used 
Bayes theorem to write the probability that the sys- 
tem is in configuration p — {p\, p2, ■ ■ ■ , pn} at time step 
p as P(pi,pi +1 ;p)P(p 1 , . . .,p/_i,p i+2 , . . .,p N \pi,pi +1 ;p); 
in which P(pi, pi + i;p) is the probability of finding sites 
I and / + 1 with energies pi and at time p, indepen- 
dently of the configuration of the remainder of the sites, 



and P{pi, . . .,pi-i, Pi+2, ■ ■ ■ , Pn\Pu Pl+x\p) is the condi- 
tional probability of finding the remainder of the sites in 
configuration {px, . . . , pi-i,pi+2 ■ ■ ■ , Pn} at time p, pro- 
vided that sites I and / + 1 have energies pi and pi+i, 
respectively. 

First, let us calculate the integral over 
{pi,...,pi-i,pi +2 ...,pN} in Eq. dAgt. From the 



definition of fl p (L), Eq. (2.4) 
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'1-2 



Q(L) 



^2f(pi' + pi'+i) + f(pi-i + pi+2) + f(pi'+pi'+i) 



U'=l 



l'=l+2 



1 

z 



i+1 



f{pv + pi'+i) - /Oi-i + pi+2) 



(A3) 



Note that, in Eq. (A3), we write f2(L) and pi, omitting time dependence appears in the probability distributions. 



the subindex corresponding to the time step p. This is so 
because we will insert Eq. ( |A3[ ) into ( |A2| , where the in- 
tegration variables are dummy, {pi, p2, ■ ■ ■ , Pn}, and the 



In fact, we have already used this notation in Eq. (A2|, 
for flp(L) but also for the functions g and /. Equation 
( A3 ) implies that 



1+1 



Y f(Pi' + Pi'+i) - f(Pi-i + Pw) 



x>=i-i 



(A4) 



where 0*(L — 2) corresponds to a system with size L — 2, without the sites / and I + 1. Equation ( A4| suggests that, 
in the large system size limit, 



dp 



1 • • • 



dp, 



1-1 



dp, 



7+2 ■ ■ ■ 



dp N 9, 1 (L)P(p 1 ,...,pi„ 1 ,pi +2 ,---,PN\pi,Pi+i;p) ~ (H* 1 (L-2)) 



(A5) 



where r p = (fi p (L)} for i — > 00, as defined in Eq. average value is calculated at time p, with the probabil- 
(3.11b). The subindex p reappears in r p , because the ity distribution P(p±, P2, ■ ■ ■ , Pn\p)- Thus, substituting 

Eq. (|A5| into (|A2l, 



)f(fii, P + Pl+x,p) 

n p (L) 



t p (g(pi, P , Pi+i, P )f(pi,p + pi+i, P )) 



(A6) 



This result will be repeatedly used in the following sections to simplify the calculations. In order to further advance, 
we have to calculate the average 



(g(pi, P , Pi+i,p)f(pi, P + Pi+i,p)} 



dpi / dpi + xg(pi,pi +1 )f(pi+ pi +1 )P{pi,pi +1 ;p), 



(A7) 



for which the distribution function P(pi, pi +1 ;p) must be known. At this point, and so as to calculate these averages, 
we introduce the local equilibrium approximation, i.e. we assume that the probability distribution P(pi, pi+\;p) inside 
the integral in Eq. ( A7 ) can be substituted by 



Ple{pi,Pi+i;p) = 



(pl, P ){pl+i,p) 



x exp 



Pi 



Pi+i 



, -(A8) 

{Pl,p) W+i,p) J 

This local equilibrium hypothesis allows us to calculate 
expressions of the averages in terms of local energy den- 
sities (pi t p) and Furthermore, it will be assumed 
(when necessary) that (p;+i, p — pi, P ) = C(L _1 ), consis- 
tently with the spatial continuum limit introduced in sec. 
[Ilj In this way, the average in Eq. ( |A7[ ) becomes a func- 
tional of the energy density field p av {x, t) , 

/ v _ / v , j-i (Pl+i,P - Pl,p) 
\Pl+l,p) - \Pl,p) + L> 



Ax 



p av (x,t) + L 1 d x p av (x,t). 



(A9) 



since 



(pi+i, P - Pi, P ) 
Ax 



d x p av (x,t) 



(A10) 



in the spatial and time continuum limits we are consid- 
ering throughout the paper. 



1. Derivation of equation ( |3.14[ ) 



Our starting point is Eq. (3.7) for the current, 



Ul, P ) = 



" / (Pl,p - Pl+i, P )f(£l,p) 



2L 
a 

2L 



n p (L) 

Tp{(Pl,p- Pl+l,p)f(£l,p)), (All) 



16 



where we have made use of Eq. ( A6 ). Therefore, we have 
to evaluate the average ((pi, p — pi+i tP )f(Hi tP )) in the local 
equilibrium approximation. Up to first order in L , 



((Pi,p ~ Pl+i,p)f(pl, P )) 



LE 



-i d x p a 



p 3 

rav 



poo poo / 

J dpi J dp i+1 (pi - pi+i) ( 1 



Pl + l \ ft A \ IP 1 +PI+ 1 



(A12) 

where we have taken into account Eqs.(A8) and (A9). By interchanging the dummy integration variables pi and 



Pi+i, an expression symmetric to the previous one is obtained. Averaging the two equivalent expressions, the above 
equation is transformed into 



((pi, P ~ Pi+i, P )f(Ei,p)) 



LE 



-L 



-xd x p 



2pf 



dpi / dpi+i(pi - pi+i) f(pi + pi+i) exp 

av "'O JO 



Pi + pl+l 



= -L 



-id x p a 



dx 
Jo 



dy(x-y) 2 f(p av (x + y))e-^ + y\ 



(A13) 



where we have made use of (p;+i iP ) ~ (pi,p) Pav( x it), 



as given by Eq. (A9|, in the large system size limit. 



Besides, in order to obtain the last equality in Eq. ( A13 1, 
we have introduced the change of variables pi = p av x, 
pi+i = p av y- We can further simplify this expression by 
going to polar coordinates (r, </>), and defining x = a 2 , 
y = 6 2 , with a = r cos </>, b = rsinc/i. Thus, the integral 
over the angle variable cp can be carried out, with the 
result 



((pl,p-Pl+l,p)f(El,p)) 



LE 



d x pa* 

3L 



drr 7 f(p av r 2 )e~ 



(A14) 

The equation above, together with Eqs. (All I, is equiv- 
alent to Eqs. (13. 141)- (13. 15b. 



2. Derivation of equation (3.17| 



Now, we will start from Eq. (3.8) for the dissipation, 



(di, P ) 



1-a 
L 

1 - a 



(A15) 

where we have used again Eq. ( |A6[ ). Similarly to what 
has been done in the previous subsection, we calculate 
now the average (pi }P [/(E/ jP ) + /(E;_i iP )]) in the local 
equilibrium approximation. In fact, we have a simpler 
case, because this average is of the order of unity in the 
large system size limit. It is easily shown that 

(Pl,p [f(El,p) + /(E;_i, p )])lE = Pav 

,h I dy(x + y)f(p av (x + y))e-^+y\ 

(A16) 



As already done when writing Eq. (A13), we have taken 
into account (pi+i, p ) ~ (pi, p ) -> p av (x,t), Eq. (|A9|, and 



introduced the change of variables pi = p av x, pi + i = 
Pnvy- Now, by introducing the same change of variables 
to polar coordinates as in the previous subsection, and 
making the integral over the angle variable, we arrive at 



(p hp [/(E iiP ) + /(E/_i iP )])le = 2 Pa 



< I drr b f(fi av r' l )e- 
h 



(A17) 



that implies Eq. (3.17), taking into account Eqs. (A15) 



3. Derivation of Eq. (4.10| 



Now, we evaluate the average in Eq. ( |4.9[ ), 

[pf+t, P + Pl,p - Pi,pPi+i,p) 



a 

T 



n p {L) ' 

■ ((pf+i,p + Pip ~ Pi, P Pi+i,p) f(£l,p)) 



(A18) 

along the same lines of the previous subsections. By in- 
troducing the local equilibrium approximation, 



{{Pi+hp + Pl, P ~ Pl,pPl+l, P ) /( s «,p)) 



LE 



:/ dx dy{x 2 +y 2 - xy)f{p av (x + y))e {x+v \ 
lo Jo 

(A19) 

where we have used again the same simplifications as 
above. After going to polar coordinates and evaluating 
the angle integral, one arrives at 



((P/+i, P + P\p ~ Pl,pPl+i, P ) /( s i, P )) 



LE 



x / drr 7 f(p av r 2 )e- r . (A20) 
Jo 
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Inserting Eq. (A20) into Eq. (A18| leads to the desired 
result p~Io|). 



Appendix B: Gaussian character of the current noise 



In the large system size limit L> 1, the current noise 



introduced in section IV B is white, 



(COM)) = 0, *')) = T o{p)5(x-x')8(t-t'), 

(Bl) 



L 



with the mobility a(x,i) given by Eq. (4.151, indepen- 



dent of L. We can introduce a new noise field £(x,i) 
by 



(B2) 



and £(x, t) remains finite in the large system size limit as 
L — s- oo, 

(^x, *)) = 0, {£(x, t)Z(x', t')) = a(p)S(x - x')S(t - t'). 

(B3) 

In this appendix, we show that all the higher-order cu- 
mulants of t) vanish in the thermodynamic limit as 
L — > oo. Let us consider a cumulant of order n of the 
microscopic noise that is equal to the n-th order mo- 
ment of £ plus a sum of nonlinear products of lower order 
moments of £. A calculation analogous to the one carried 
out for the correlation {£i, P £i>, P <) shows that the leading 
behavior of any moment is of the order of L , which is 
obtained when all the times are the same. Therefore, the 
moment (ji, p ji>, p > • • ■ Jj(n) »c»o) gives the leading behavior, 
of the order of L" 1 for p = p' = ■ ■ ■ = p( n '; any other 
contribution to the cumulant is at least of the order of 
L- 2 . Therefore, 



{jl,pjl',p> ' ' ' j;("),p<")) ~ L 1 Tp(Cl p }5ll>5l>,l>> ■ ■ ■ 5j(n-l) ) J(n)flp ) p'5p/,p/' • • • <5p(»-l) ip (n) , 



(B4) 



where (Ci p ) is a certain average which remains finite in the large system size limit a s L — > oo. In the continuum limit, 
each current introduces a factor L 2 /t p due to the scaling introduced in Eq. (3.11a), and 



-2n-l 1-n 



(C(x,t)) 



S(x - x')S(x' - x") ■ ■ ■ 8{x^ n -^ - x^) 



rfH(t - t')6(t' - t") ■ ■ ■ S^"" 1 ) - tW) 

L3(n-1) 



(B5) 



i.e. 



t)i{x\ t 1 )--- ttx [n \t [n) )) = L 3 - 2n (C(x, t))5(x - x')S(x' - x") ■ ■ • - x (n) ) 

xS(t - t')S(t' - t") ■ ■ ■ 5(& 1 -^ - iW), 



(B6) 



where (C(x, t)) is the (finite) leading behavior of (Ci :P ) in the large system size limit L 1, under the local equilibrium 
approximation int roduc ed in section |IV A| We have also taken into account the relationship between Kronecker and 
Dirac deltas, Eq. (4.16|. Going to the rescaled, of the order of unity, noise £, 



t)i(x', If)--- i (n) )) ~ L 3 ( x -f ) (C(x, t))8{x - x')5{x' - x") ■ ■ ■ 5{x {n -^ - x [n) ) 

xS(t - t')8(t' - t") ■ ■ ■ 5(t (n - r > - iW). 



(B7) 



Therefore, in the limit as L — > oo, 
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which completes the proof. A similar calculation can be 
carried out for the dissipation noise, that is also gaussian 
in the thermodynamic limit. Nevertheless, it will not 
be presented here, since it is subdominant as compared 
to the current noise for large system sizes L ^> 1, and 
therefore has been neglected. 
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